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Abstract 

Isotropic models in loop quantum cosmology allow explicit calculations, thanks 
largely to a completely known volume spectrum, which is exploited in order to write 
down the evolution equation in a discrete internal time. Because of genuinely quan- 
tum geometrical effects the classical singularity is absent in those models in the 
sense that the evolution does not break down there, contrary to the classical situ- 
ation where space-time is inextendible. This effect is generic and does not depend 
on matter violating energy conditions, but it does depend on the factor ordering of 
the Hamiltonian constraint. Furthermore, it is shown that loop quantum cosmology 
reproduces standard quantum cosmology and hence (e.g., via WKB approximation) 
to classical behavior in the large volume regime where the discreteness of space is 
insignificant. Finally, an explicit solution to the Euclidean vacuum constraint is dis- 
cussed which is the unique solution with semiclassical behavior representing quantum 
Euclidean space. 



1 Introduction 



General relativity is a very successful theory for the gravitational field which is well tested 
in the weak field regime. However, it also implies the well-known singularity theorems QXJ] 
according to which singularities and therefore a breakdown of this theory are unavoidable 
provided that matter behaves in a classically reasonable manner (i.e., fulfills energy con- 
ditions). In fact, observations of the cosmic microwave background demonstrate that the 
universe was in a very dense state a long time ago, which classically can be understood 
only in models which have an initial singularity. 

However, if the universe is small and the gravitational field strong, the classical de- 
scription is supposed to break down and to be replaced by a quantum theory of gravity. In 
early attempts it was proposed that the singularity could be avoided by coupling classical 
or quantum matter fields which violate the strong energy condition and thereby evade the 
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singularity theorems 0, |3|, £|, [^]. Another approach modifies the field equations by adding 
higher curvature terms interpreted as the leading order corrections of quantum gravity 
H [7], |[]. However, the first approach uses matter effects rather than those of quantum 
gravity, and in the second, more and more corrections are needed the closer one comes to 
the classical singularity. Furthermore, the procedure of truncating the series of perturba- 
tive corrections (involving higher derivatives), and treating all solutions to the resulting 
equations of high order on the same footing as those of the lowest order ones, is inconsistent 
||. Eventually one still needs a full non-perturbative quantum theory of gravity in order 
to understand the fate of the singularity. 

Lacking a full quantization of general relativity, an early approach to quantum cosmol- 
ogy was to perform a symmetry reduction to homogeneous or isotropic models with a finite 



number of degrees of freedom and to quantize afterwards |T0|, O] ; this will be called "stan- 



dard quantum cosmology" in the following. But in general, those models do not avoid the 
classical singularity (even the meaning of this phrase is not clearly understood: an early 
idea was to impose the condition that the wave function ip(a) vanishes at the singularity 
a = PH |, but this is insufficient with a continuous spectrum of the scale factor a JI2]] 
which is always the case in standard quantum cosmology), and more severely quantum 
mechanical methods are used for the quantization which are not believed to be applicable 
to a full quantization of general relativity. Therefore, it is not clear to what extent the 
results are relevant for quantum gravity. 

In the meantime, a candidate for quantum gravity has emerged which is now called 
quantum geometry [ 13| , |14| . A key success of this approach is the derivation of a discrete 
structure of space which is implied by the discreteness of spectra of geometric operators 
like area and volume [15], [16|, [17j]. It also leads to a new approach to quantum cosmology 
T8|] : By reducing this kinematical quantum field theory to homogeneous or isotropic states 



using the general framework for such a symmetry reduction at the quantum level [|T^] one 
arrives at loop quantum cosmology |20[| . Due to this derivation loop quantum cosmology 
is very different from standard quantum cosmology; e.g., the volume is discrete which 
is inherited from quantum geometry. As a consequence the Hamiltonian constraint, which 
governs the dynamics, of those models can be written as a difference equation rather 
than the differential, Schrodinger- or Klein-Gordon-like evolution equation of standard 
quantum cosmology; thus, also time is discrete P5| . One may expect that the discrete 
structure of space-time, which is most important at small volume, will have dramatic 
consequences for the appearance of a singularity^. As a loose analogy one may look at the 
hydrogen atom: classically it has a continuous family of orbits leading to its instability, 
which is quantum mechanically cured by allowing only a discrete set of states. Of course, 
one also has to determine radiation loss or transition rates to judge stability, i.e. one has 
to take into account the fully dynamical situation. 

In fact techniques developed for a quantization of the full theory [f24]] can be used to show 
in isotropic models that the inverse scale factor, whose divergence signals the singularity in 
the classical model, can be quantized to a bounded operator with finite eigenvalues even in 
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states which are annihilated by the volume operator This is already a first indication 
that close to the classical singularity the system is at least better behaved in a quantum 
theory of geometry, which is to be confirmed by studying the dynamics of these models 
as governed by the Hamiltonian constraint equation. The results of this analysis have 
been communicated in |26[] and are described in detail in the present paper: whereas the 
classical singularity presents a boundary of space-time which can be reached in finite proper 
time but beyond which an extension of space-time is impossible, the quantum evolution 
equation does not break down at the classical singularity provided only that we choose the 
appropriate factor ordering of the constraint. The conclusion is, however, independent of 
the particular matter coupled to the model, and so does not rely on special forms of matter 
violating energy conditions. We will describe these conclusions in detail in Section (| after 
reviewing the formalism of isotropic loop quantum cosmology (Section |2|) and computing 
the action of the Hamiltonian constraint for spatially flat models (Section ^ the constraint 
for positively curved models can be found in App. IB]). After discussing the fate of the 



singularity in Section [4.2| we will see in Section ^]3| that we simultaneously have the correct 
semiclassical behavior at large volume, and present in Section |5.1| an explicit solution to the 
vacuum Euclidean constraint which gives quantum Euclidean space. A comparison with 
other approaches and ideas put forward to avoid or resolve the singularity can be found 
in Section |6| before our conclusions where we will also discuss in which sense, compared to 
the classical singularity theorems, the singularity is absent in loop quantum cosmology. 



2 Isotropic Loop Quantum Cosmology 

A calculus for isotropic models of loop quantum cosmology can be derived from that of 
quantum geometry along the lines of the general framework for a symmetry reduction 



of diffeomorphism invariant quantum field theories |]19| . This allows us to perform the 
symmetry reduction at the quantum level by selecting symmetric states. In the connection 
representation, those states are by definition only supported on connections being invariant 
with respect to the given action of the symmetry group. They can be embedded in the full 
kinematical Hilbert space as distributional states. 



Specializing this framework to isotropy |2(J leads to isotropic states supported on 



isotropic connections of the form A % a = 0}co^ = cA)u^. Here, the internal S'f/(2)-dreibein 
A/ = A}Tj is purely gauge {jj = —\oj are generators of SU{2) with the Pauli matrices <x,) 
and lo are left-invariant one-forms on the "translational" part N (isomorphic to M. 3 for the 
spatially flat model or SU(2) for the spatially positively curved model) of the symmetry 
group S = N x SO (3) acting on the space manifold E. Orthonormality relations for the 
internal dreibein are 



A} A/ = 5j and e ijk A)A J jA k K = e IJK . (1) 

For homogeneous models, the nine parameters (ft) are arbitrary, and for isotropic models 
c is the only gauge-invariant parameter. A co-triad can be expressed as e\ = a\io l a = 
aA)uJa with the scale factor \a\ (in a triad formulation it is possible to use a variable a 
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which can take both signs even though the two corresponding sectors are disconnected 
in a metric formulation). Using left-invariant vector fields Xj fulfilling u/(Xj) = 8j, 
momenta canonically conjugate to A l a are densitized triads of the form Ef = pfXf = pA^Xf 
where p = sgn(a)a 2 . Besides gauge degrees of freedom, there are only the two canonically 
conjugate variables c and p which have the physical meaning of curvature and square of 
radius, and are coordinates of a phase space with symplectic structure 

{c,p} = hn (2) 



(k = 8ttG is the gravitational constant and 7 the Barbero-Immirzi parameter; see |25[ for 
an explanation of the correct factor | which is missing in J2"T|). 

Any gauge invariant isotropic state, being supported only on isotropic connections, can 
be expressed as a function of c. An orthonormal basis of such functions is given by the 
usual characters on SU (2) 

sin (7 + i)c . , , 



sm |c 



together with C-i — (v2sin|c) 1 and 



cos(? + ~)c , , 

= 1 2 > ie|N . (4) 



sm 



Gauge non-invariant functions are given by A)xj and A^Q where A} is the internal dreibein 
providing pure gauge degrees of freedom. 



The states Xji Cj are a ls° eigenstates of the volume operator [^] with eigenvalues 



v 3 = {iii)^yu + \)u + i)- (5) 

Because j can take the value — | as label of Q the eigenvalue zero is threefold degenerate, 
whereas all other eigenvalues are positive and twice degenerate. An extension of the volume 
operator to gauge non- invariant states is done by using the relation [A}, V"] = which 
follows from gauge invariance of the volume . 



Because the basic multiplication operator is the point holonomy hi := exp(cA}rj) = 
cos(|c) + 2 sin(|c)A}rj we also need the action of cos |c and sin |c. This can be obtained 
in the connection representation (^|), (|j) by using trigonometric relations, but it is easier 
first to introduce a new orthonormal basis of the states 

. . exp(iinc) , , 

|n>:= ^ 2 r - ^ Z ( 6 ) 
V2 sm 

which are decomposed in the previous states by 

\n) = 2-3 (Ci(H-i) +«sgn(n)xi ( |„|_i)) 
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for n ^ and |0) = The label n, which will appear as internal time label below, is 



the eigenvalue of the dreibein operator [21 



P 



I^p (4- + |cot|c 



V dc 



On these states the action of cos |c and sin |c is 



cos(~c)|n) = | (exp(|zc) + exp(-|zc)) 



exp(pnc) 

y2 sin(|c 

l. 



i(| n + l) + („_!)) 



(7) 



sin(ic)|n) = — |i (exp(|ic) — exp(— \icj) 2 - — — \i{\n + 1) — \n — 1)) (8 



a/2 sin^c) 



and that of the volume operator is 



V\n) = Vu lnhl) \n) = (h*p)MM-l)M(M + l)|n> 



(9) 



Together with the volume operator the trigonometric operators establish a complete cal- 
culus for isotropic cosmological models and more complicated operators can be constructed 
out of them. As a first application the expression mjj := qjj /-y/det q = |a| _1 5/j of the in- 
verse scale factor has been quantized [^5j using techniques developed in quantum geometry 



in order to quantize co-triad components, which are not fundamental variables ||24j| . The 
same method, which also leads to densely defined quantizations of matter Hamiltonians 
27] , results in a bounded operator quantizing the inverse scale factor 



mu 



16(7*1 



cos^cj V V cos(|c) 



sin(|c) V V sin(^c) 



(10) 



— 5jj ( sin(^c) V V cos(^c) — cos(|c) V V sin(^c] 



which has been studied in detail in 25 . We will later need the action 



sin(|c)Kcos(|c) - cos(|c)Vsin(|c)) \n) = \i (v"i ( | n+lhl) - Vi {]n _ lhl) J \n) 



|2Sgn(nj 



V iH - Vj w _i) |n> (11) 



(for n = the value of V-\ is understood to be zero). Important for the results of the 
present article is that the state |0) — C-i is annihilated by both the volume operator and 
the inverse scale factor. This would, of course, be impossible in a classical theory and is a 
purely quantum geometrical effect. (Since the singularity a = is not part of the classical 
phase space, one has to extend the inverse scale factor appropriately which is done here 
formally by sgn(a) 2 /a [25|.) Our later considerations crucially depend on the fact that 
all metrical operators, and therefore all matter Hamiltonians which in some way always 
contain metric components, annihilate the state |0). 
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3 Hamiltonian Constraint for Isotropic Models 



For homogeneous models |20] the Euclidean part (which is the full constraint in Euclidean 
signature if 7 = 1) of the Hamiltonian constraint is given by (the lapse function is irrelevant 
for cosmological models and set to be one) 

H£L = -K-Ueiia^e^FjjE^Ei 

= n- 1 det(a})- 1 (e ufc cf J 0^ J - 4<t>jPjPk + Mjv]vl) (12) 

[F\j are the curvature components of the connection A\ and cfj are the structure constants 
of the symmetry group). Here, we had to choose a relative sign for the constraint in the two 
different orientations of the triad. Classically, both orientations are disconnected because 
one has to require a non-degenerate triad. However, this is no longer necessary in a quantum 
theory, and we will in fact see that an evolution through degenerate metrics is possible. 
Therefore we have to choose the relative sign which we did by using the determinant 
of the co-triad instead of the metric, which would always be positive. The quantization 
techniques of |24| directly apply only to the convention in fll2|) . We will see that one can 
transform between the two choices after quantization, however not unambiguously due to 



special features at the classical singularity; see the remarks following (|22|) . 



3.1 The Classical Constraint 

In the present context of isotropic models we are interested only in special homogeneous 
models which can be further reduced to isotropy. These are the two Bianchi class A models 
given by the structure constants cfj = (Bianchi type I) which lead to the isotropic 
spatially flat model or cfj = efj (Bianchi type IX) leading to the isotropic model with 
positive spatial curvature. The third isotropic model, which has negative spatial curvature, 
can only be derived from a class B model and so is not accessible in the present framework. 
Inserting isotropic connection and triad components into fll2|) yields 



#( E ) = Q K - 1 a- 3 (2T - c)cp 2 = 6k- 1 (2T - c)csgn(p)^/\p~\ (13) 

where T = for the flat model and r = ~ for the positively curved model. 

The parameter T also determines the spin connection compatible with a given triad, 



which is given by [28 



For homogeneous triads e" = a\X^ and inverse co-triads e l a = a}uj^ this specializes to 
which using the Maurer-Cartan relations d[ a u>^ = —^Cj K uj£(j£ yields 

■pi 1 Ajk/„K _ 1 1„_L J ' K J, „l \, ,1 -pi , ,1 

1 a ~ 2 C ( C U a j a K + 2 C KJ a j a k a I a L) UJ a — 1 I^a ■ 
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In isotropic models the co-triad has the special form a) = aA] which implies a\ = a X A\ 
and leads to 

O = \<$AM + l Aje JKM ^iL 

using the ortho normality relations ([!]) for the internal dreibein A}. Now we use that 
for Bianchi class A models the structure constants have the form cfj = euLn LK where 

With these structure constants we finally arrive at 



/; . <> = n^ K '5 LK is a diagonal matrix with n 1 = for Bianchi I and n 1 = 1 for Bianchi IX. 



= \(n l +n 2 + n 3 - 2n {I) )A i I = TA) (14) 

where the constant T which specifies the isotropic model has been defined above. In 
particular, we see that the spin connection vanishes for the spatially flat model, and so the 
Ashtekar connection A l a = T l a + ^K^ is proportional to the extrinsic curvature, whereas in 
the positively curved model it has an extra term given by the intrinsic curvature of space. 

Knowing the spin connection we can express the integrated trace of the extrinsic cur- 
vature in terms of the isotropic variables (c,p): 

K := J d 3 xKiE« = 7" 1 J d 3 x(^ - VjEf 

= 7" 1 (0/ - typl = 37^(c - I> =: Kjpl . (15) 
We also define the isotropic extrinsic curvature component k by kA) := K\ = 7~ 1 (c — T)A} 



such that K = 3kp. As proposed in [24], we will use this quantity in order to quantize the 



Hamiltonian constraint for Lorentzian signature exploiting the relation 

K = ! 7 ~ 2 {# (E) ,^}sgn(det4) 



which can easily be verified here for the isotropic Euclidean part (O) of the constraint with 



the symplectic structure (0). This can then be inserted into the Lorentzian constraint 

H = -H {Fl) + P (16) 

with 

P := -2(l + 7 2 ) K - 1 det(4)- 1 Kf / ^ ]J E/E/ 

= -(1 + 'r, :i (A-;/^///,/ - KiK'pipj) 

= -6(l + 7 2 )fiT 1 a- 3 A;y = -6(l + 7- 2 )/?- 1 (c-r) 2 sgn(p) v 4pj (17) 
to yield the Lorentzian constraint for isotropic models. 

3.2 Quantization 



According to |22[, the Hamiltonian constraint for homogeneous models can be quantized 
along the lines of the full theory |M] if the special requirements of the symmetry are taken 
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into account. One arrives at the Euclidean part (note that this corresponds to the relative 
sign for the two different triad orientations as chosen above; see the discussion following 
equation (|T2|)) 



ff( E )=4i( T ^)- 1 ^e / ^tr(Mj/i7V^j]^[^ 1 .^]) 



(18) 



UK 



where the holonomy operator h[i > j] depends on the symmetry type and is defined by 



IK, 



K 



In contrast to we quantized the Poisson bracket of (f)} and the volume to —hK[h K , V"] 
which is completely along the lines of the full theory. Although in homogeneous models it 



is possible to use the simpler expression [/i^V], which has been done in [22|, this is not 



advisable as can be seen from a quantization of the inverse scale factor | 25| . Therefore, we 
use here the quantization which is closer to that in the full theory. 

Using the extrinsic curvature, the Lorentzian constraint operator can be written as 



H 



.#(E) + p 



_£CB) _ 8i(1 + 7 2 )«- 1 ( T ^)-3 e IJK tr fh^hj 1 , k]hj[h-\ K]h K [h- K \ U])(19) 



UK 



with (up to ordering ambiguities which will be discussed below) 

K = -\ %1 - 2 h- Y [# (E) , v] sgn(deTa}) . 



(20) 



Since K appears in P within a commutator with holonomies (corresponding to Poisson 
brackets in the classical expression), the sign sgn(deta)) in ( p0|) is important even though 
P is quadratic in K. 

From the homogeneous operators we can always derive the isotropic ones by inserting 
holonomies hi = cos(|c) + 2sin(|c)A}rj and the isotropic volume operator. Using the 
dreibein relations ([]]) for A} and tr(Aj) = one can then take the trace in order to arrive 
at an operator composed of cos(|c), sin(|c) and V. For the Hamiltonian constraint we 
need 

hjhjhfhf = cos 4 (|c) + 2(1 + 2e /J i %0 sin 2 (|c) cos 2 (|c) + {25 u - 1) sin 4 (i C ) 
+4(A 7 - Aj)(l + 8u) sin 3 (i C ) cos(ic) 

which has been computed using A/Aj = ^€u k Ak — \8ij. Similarly, we have 

hxih^jV] = V — cos(t;c)V cos( |c) — sin( \c)V sin( tjc) 

— 2A^ ^sin( |c) V cos ( ~ c) — cos(tjC)U sin( |c) J . 



8 



These formulae are sufficient in order to derive the Euclidean part of the Hamiltonian 
constraint for the spatially flat isotropic model 

#( E ) = 9Qi(-fKll)- 1 sin 2 (±c) cos 2 (±c) (sin(±c)V cos(fc) - cos(±c)y sin(Jc)) (21) 

with action 

H^\n) = 3( 7 ^)- 1 sgn(n) (Vi |n| - V^) (|n + 4) - 2\n) + \n - 4)) (22) 

using (11). Here we see that one can transform to the other sign convention in (12) 
simply by dropping sgn(n). At this place it looks unambiguous since for n = we have 
Vi| n i — V r i| n i 1 = 0, but note that the splitting into the sign and the difference of volume 
eigenvalues in (|TT|) is ambiguous (since we use sgn(0) := 0, its prefactor is not uniquely 
defined; in (|TTD it has just been extended from the general expression for positive and 
negative n). In quantum geometry only the constraint with sign convention as in (|T2| ) can 
be quantized directly and appears much more natural. 

Next we can build the extrinsic curvature operator using the Euclidean part of the 
Hamiltonian constraint. Since we chose a non-symmetric ordering for the Euclidean con- 
straint, we obtain from (EH) a non-symmetric extrinsic curvature operator 



K\n) = I^V^h-Vii^! 



Vi ( i n+4 |_ i) - Vu M _ n ) |n + 4}- (Vifu-n -Vu\ n M X) ) |n-4) 



X 1 "i(|n+4|-l) - K |(N-l)y l'*-r-*/ - y±(\n\-l) ~ "i(| rl -4|-l) / 

= : §i£(/£|n + 4)-£+|n-4>) (23) 
where we defined the coefficients 

K± : = Tl2(7^)~ 3 (V hln] - V, lnhl ) (v UlnT , hl) - V hQnhl) ) (24) 

which fulfill JCZ n = — and are approximately given by /Cj ~ n for large \n\. It is 
possible to choose a symmetric ordering of fC without changing the original ordering of 
H( E \ but this is not necessary since we are interested here only in the constraint, which 
need not be symmetric (in fact, the Euclidean part must not be symmetric as we will see 
below). Nevertheless, we will see shortly that not all orderings for 

K\n) := a/C+(l-a)/C t = \i$(K-\n+4)-K+\n-4)) , K± = a/C± + (l-a)/Cj T4 (25) 

with ordering parameter a G M. are allowed. Using K in some given ordering we obtain 
the potential term of the Lorentzian constraint 

P = -8*(l + 7 2 K%/|r 3 ^e /J *tr(M^ 

UK 



x (sin(-^c)K cos(|c) — cos(^c)K sm(^c)^ (^sm(^c)V cos(^c) — cos(^c)V sin(^ 



2 C > 



9 



We already commented on the ordering of K, and it will be seen to be crucial to order the 
operator containing V, which quantizes the triad components, to the right. 
For P we need the action 

(sin(±c)iTcos(±c) - cos(§c)ifsin(±c)) |n) =: -\ll{k~\n + 4) - k+\n - 4)) (26) 
introducing 

*£ = &K&i-Kti) (27) 

with fc+ = kZ n . Inserting fl2"E|), fl2"T| ) and (^) into the action of P demonstrates that 
expressions for the action of the isotropic Lorentzian constraint can be quite cumbersome, 
but can be computed explicitly thanks to the completely known volume spectrum (the 
equation simplifies at large volume where the coefficients k^ are approximately one, as 
used in Section |5.2|) . Fortunately, for the later discussion we will not need the explicit 
action but only the crucial fact that we can order K such that the coefficients k^ are 
non-vanishing for all n. In our original operator (^3|) the coefficients /C^ and JC± 2 vanish 
leading to k^ = which will be seen in the next section to give a singular evolution. One 
can easily remedy this by ordering K symmetrically which amounts to replacing /Cj by 
\{fCn + K-n^) an d results in coefficients k^ which never vanish. From now on, we use the 
ordering 

such that 

K n = 2 Y^n + ^n=F4) • 

Finally, we can compute the action of P by using the previously derived operators: 

P\n) = |(l + 7 - 2 )(7^)- 1 sgnH(^ w -^ |nhl ) (28) 

x {KK+M + 8) - (K K+i + KK-i)\ n ) + KK-i\ n - 8 )) • 
4 Evolution 

Now we have all ingredients for the explicit form of the Hamiltonian constraint equation 
and can discuss the evolution it governs. For simplicity, we will write down the following 
formulae for the spatially flat model and will only comment on possible differences in the 
spatially positively curved model, but most qualitative results apply to both cases. 

4.1 Discrete Time 

We will use the triad coefficient p = sgn(a)a 2 as internal time which classically makes sense 
only for positive p. In particular, the evolution breaks down at the classical singularity 
p = so that the two branches p > and p < are disconnected. (Note that p and 
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— p result in the same metric?], but are not identified by a gauge transformation since the 
gauge group is £0(3) rather than 0(3). Factoring out the large gauge transformation 
p — ► — p is not allowed since, in particular, it results in a conical singularity at p = in 
the extended phase space which is used in the quantum theory.) This setting enables us to 
discuss the fate of the singularity in quantum cosmology by studying the evolution close 
to p — 0. For this we need to write the constraint equation H\s) = for a history \s) 
as an evolution equation which can only be done in a dreibein, rather than connection, 
representation since we are using a metrical expression as internal time [29, 13]. According 
to |23H , a dreibein representation is defined by expanding \s) = J2n s n\ n ) an d using the 
coefficients s n as a wave function in the dreibein representation. Because of the discreteness 
of geometric spectra a state in the dreibein representation is a function on a discrete set 
given by Z here, and using a metrical internal time implies a discrete time evolution p3[ . 
With our explicit expressions fl22|) , ( p8]) for the Hamiltonian constraint we can write down 
the difference equations governing the evolution of isotropic models. 

In the dreibein representation, the Euclidean part of the constraint acts as 

(H {E) \s)) n = 3( 7 ^)- 1 [sgn(n + 4)(^ | „ +4| -Vi |n+4hl ) Sn+ 4 

- 2sgn(n)(Vi w - Vu nhl )s n + sgn(n - 4)(Vi| n _ 4 | - Vi | n _ 4 |-iK-4 



and P as 



2\-l 



sgn(n + 8) ( Vi 



\n+8\ 



T \n+8\- 



|(l+ 7 -)( 7 ^) 

- sgn(n) (Vi H - Vi| nhl ) {k~k+ +i + fc+fc"_ 4 )s n 
+ sgn(n - 8) (Vi| n _ 8 | - Vi|„_ 8 |_i) K-sK-4 s n-8 



l ) k n+8 k n+i s 



n+8 



In a realistic cosmological model we also need matter which enters the evolution equa- 
tion via its matter Hamiltonian, and we may include a cosmological term. The precise 
form of the matter and its quantization is not important here, and we will build it into 
our description by using states s n (4>) in the dreibein representation which are functions 
of the matter degrees of freedom <p. Since matter and gravitational degrees of freedom 
are independent and so commute prior to imposing the constraint, we will get a matter 
Hamiltonian acting on s which is diagonal in the gravitational degree of freedom n for 
usual matter (for the rare cases of matter with curvature couplings our discussion has to 
be adapted appropriately; note that our starting point is not an effective Hamiltonian with 
possible higher curvature terms which would have to be derived only after quantization): 
H^\n) <S> |0) =: \n) <8> Hj,(n)\<f>). Important for what follows is also that all states |0) ® \<f>) 
are annihilated by the matter Hamiltonian, which is a consequence of special features of 
quantum geometry. More precisely, all terms in a matter Hamiltonian contain metric com- 
ponents in order to have a scalar density which is integrated to the Hamiltonian; and a 



2 The author thanks Y. Ma for making him aware of this issue. 



11 



quantization of those metric components along [27 leads to an operator annihilating |0) 



and so -£^(0) = 0. This is analogous to the inverse scale factor which also annihilates 
|0) and is possible only in quantum geometry because in a matter Hamiltonian both the 
metric and the inverse metric can appear which cannot simultaneously be zero classically 
(see JZS) for a detailed discussion). 

Collecting all ingredients we arrive at the evolution equation 

|(1 + 7~ 2 ) Sgn(n + 8) (Vi| n+8 | - Vi|n+8|-l) k n+8 k n+4Sn+8(<P) 

- sgn(n + 4) (Vi | n+4 | - Vi| n+4hl J s n+4 (0) 

-2 sgn(n) (Vi H - Vi| Bhl ) (|(1 + T 2 ){KK + a + KK-a) ~ l) 

- sgn(n - 4) (Vi , n _ 4 | - Vi | n _4 hl ) s n _ 4 (0) 

+K 1 + 7~ 2 ) sgn(n - 8) (Vi| n _ 8 | - Vi| n _ 8hl ) ^„ 8 A;~_ 4 s n _ 8 (0) 
= -§7/dpSgn(n)i^(n)s ri (0) (29) 



which, as anticipated in p3| , is a difference equation for the coefficients s n (0) in the discrete 
label n (our discrete time) of order 16 (in p3f the order is higher since the Poisson brackets 
with the volume have been quantized differently). Note also that due to k+ = kZ n the 
equation is symmetric under time reflection n \— > —n, provided the matter Hamiltonian 
fulfills H^n) = H^-n). 

The gravitational part of the evolution equation is quite complicated and an explicit 
solution is possible only in simple cases like the Euclidean vacuum equations discussed 
below. But the equation is amenable to a numerical analysis because it is a difference 
equation, and given some initial data we can compute subsequent components of the wave 
function (in a numerical analysis one has to be aware of possible unstable solutions P0fl). 
However, a recursive computation is possible only as long as the highest order coefficient, 
sgn(n + 8) {V\ n +8\/2 — Vj n+8 |/ 2 _i) k^ +8 k^ +i , does not vanish. As discussed earlier, we use 
an ordering of the extrinsic curvature such that k+ never vanishes, but the rest of the 
coefficient is zero if and only if n = —8. This means that, starting at negative n, we can 
determine components s n of the history s only up to n — — 1, and the coefficient for n = 
is not determined by the evolution equation. Because the volume vanishes in the state |0) 
it seems that as in the classical theory there is a singularity in isotropic models of loop 
quantum cosmology in which the evolution breaks down. This, however, is not the case as 
we will show now. 



4.2 Fate of the Singularity 

We assume that we are given enough initial data for negative n of large absolute value 
in order to specify all initial conditions for the difference equation of order 16, i.e. we 
know the wave function s n (<p) at 16 successive times n to n + 15. From these values 
we can compute all coefficients of s n (<p) for negative n, but for n = —8 the highest order 
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coefficient in vanishes. So instead of determining sq(4>) the evolution equation leads to 
a consistency condition for the initial data: We already now s n (4>) for all negative n using 
(p9|) for n < —9, and using the equation for n = —8 gives an additional condition between 
S-i6(0)) •5-12(0), s-s(0) an d s_4(</>). Such a consistency condition is not a problem because 
it serves to restrict the initial data which is welcome in order to reduce the freedom, and it 
can be used in order to derive initial conditions from the evolution equation which uniquely 
fix the semiclassical branch of a wave function However, we now lack an equation 

which would give us Sq(4>) in terms of the initial data seemingly leading to a breakdown of 
the evolution as mentioned above. But the situation is much better: We cannot determine 
so(4>) because it drops out of the evolution equation, but it also does not appear in the 
equations for n > —8. Therefore, we can compute all coefficients for positive n and in this 
sense we can evolve through the classical singularity. At this point we used the crucial fact 
H ^(0) = 0, which naturally holds in quantum geometry; otherwise, s o (0) would enter the 
evolution equation via the matter part. 

Of course, in order to determine the complete state \s) we also need to know So which, 
as we will show now, can be fixed independently of the evolution. First we note that there 
is always a trivial, degenerate solution to the constraint equation given by s n (<p) = 5o, n So(0) 
which is completely supported on degenerate metrics and which is a true eigenstate of the 
constraint with eigenvalue zero. In the vacuum case, it corresponds to the classical solution 
p = which is of no physical interest, but with matter there can be no classical analog 
because this solution arises only due to H^(0) = whereas the classical matter Hamiltonian 
usually diverges for p = 0. All other solutions are orthogonal to the degenerate solution and, 
therefore, must have So(0) = demonstrating that in an evolving solution this coefficient is 
fixed from the outset and we can determine the complete state using the evolution equation 
(in a given solution to the constraint equation there can be an arbitrary admixture of the 
degenerate state, but it does not affect the solution at non-zero n). After this discussion, 
one can absorb the sign factors in (|29| ) into the wave function, which had been done in 



31f for simplicity. This is free of ambiguities here for evolving solutions where we use 
the condition s — 0; compare the remarks after (|22|). 

We have now shown that the evolution equation (|2§| ) does not break down at the 
classical singularity, and in this sense there is no singularity in loop quantum cosmology. 
But in general a state will be supported on the degenerate states | ±1) in which the volume 
vanishes. Although this may look problematic, the inverse scale factor, whose classical 
divergence is responsible for the curvature singularity, remains finite in these states. This 
feature of quantum geometry, which also was very crucial in our proof of the absence of the 
singularity because it implied a vanishing matter Hamiltonian in |0), is the fundamental 
deviation from classical geometry leading to the consequences discussed in the present 
article. While we used this general property of the matter Hamiltonian, the precise form 
of matter is irrelevant and so our conclusion remains true for any standard type of matter 
with or without a cosmological constant, and also for the spatially positively curved model. 

On the other hand, the factor ordering of the constraint is very crucial, for in a different 
ordering the coefficient Sq(4>) would not completely drop out of the evolution and we would 
not have the degenerate solution. Ordering the triad components to the left rather than 
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to the right results in a coefficient sgn(n)(^ n |/ 2 — V\ n \/2-i) of s n+ k replacing all sgn(n + 
k)(V\ n+ k\/2 — V\n+k\/2-i) m dUD- The highest order coefficient then vanishes first for n = 
so that s$ remains undetermined and will not drop out of the equation for positive n. 
This results in si 2 , ... depending on sg, and we no longer have a solution similar to the 
completely degenerate state s n = so^no- Also, we must not use a symmetric ordering since 
in this case the highest order coefficient sgn(n)(Vj ra |/ 2 — VL|/2_i) + sgn(n + 8)(Vf n+ g|/2 — 
Vj n+8 |/2-i) vanishes for n = —4 and s 4 , s 8 , ... remain undetermined by initial data (for 
different reasons not to use a symmetric ordering of the Hamiltonian constraint in quantum 
general relativity see fl32|). Thus, a non-singular evolution of the observed kind, which is 
possible only for one of the three standard orderings (triads to the left or right, and the 
symmetric ordering) may be used as a criterion to fix the factor ordering ambiguity of the 
Hamiltonian constraint. The ordering derived here corresponds to the one chosen in |24 
for the constraint in the full theory. 

There are non-symmetric orderings of the constraint for which the highest (and lowest) 
order coefficient never vanishes. With such an ordering there would be no state correspond- 
ing to the classical solution p = in the vacuum case, and a general evolving solution would 
be supported on the degenerate state |0) with the s n depending on s . Since this state 
plays a special role even kinematically (recall that quantizations of both the scale factor 
and its inverse annihilate it [p5||), a special behavior of any evolving state, like the or- 



thogonality to it described above, is preferable. More importantly, a vanishing highest 
order coefficient implies a consistency condition which poses initial conditions on evolving 
states. As we will see in Sec. [5J], the unique state corresponding to flat Euclidean space 
only results with this condition (see also [[JI[] for the case with matter). Non- vanishing 
highest order coefficients are also obtained for a symmetric ordering if one first transforms 
to the alternative sign convention in (|22[). But recall that this transformation is not free of 
ambiguities right at the value n = which is important for a discussion of the singularity. 
Thus the choice of this ordering is problematic. The scenario of [^, ^TJ and the present 
paper is realized in only one ordering of the constraint, with triads to the right, which will 
always be used from now on. This ordering has been used previously in order to derive 
consistency of the formal constraints P5[. (However, the Chern-Simons state, which has 



been found as a solution to the formal Euclidean Hamiltonian constraint with a positive 
cosmological constant, requires the opposite ordering P3[ . But since there is no physical 
correspondence of this state it is not necessary to find it as a solution, and so its disap- 
pearance does not present an argument against an ordering; for the Lorentzian constraint 
it disappears, anyway.) 

Intuitively, we have the following picture of an evolving universe: For negative times n 
of large absolute value we start from a classical universe with large volume. It contracts 
(V(| n |_i)/2 decreases with increasing negative n) to reach a degenerate state of zero volume, 
classically seen as a singularity, in which it bounces off in order to enter an expanding 
branch and to reach again a classical regime with large volume. The change of sign in p 
during the bounce means that the universe "turns its inside out" there p5(| . A possible 
recollapse and an iteration of this behavior depends on the matter content, but for the 
evolution close to the singularity matter is irrelevant. What remains to show is that for 
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large volume we have in fact the correct semiclassical behavior, to which we turn now. 
4.3 Semiclassical Regime 

Our evolution equation (|29| ) is of the order 16 for the spatially flat model and even 20 for 
the positively curved model, and so there are many independent solutions (in fact, there 
are infinitely many independent solutions if we take into account the matter degrees of 
freedom, but we are mostly interested in the freedom coming from the sole gravitational 
variable n). Of course, a classical cosmological model does not have so many independent 
solutions, and so most of the quantum solutions cannot have a classical counterpart. In the 
present section we investigate the conditions for a solution to have a semiclassical branch. 
(This issue is discussed in more detail in |]3T|j .) 

We first have to define conditions for a semiclassical regime. Obviously, the volume 
should be large compared to the Planck scale and components of the curvature should be 
small. Moreover, continuous space-time has to be a very good approximation to the discrete 
space and time of quantum geometry. The first condition of large volume is straightfor- 
wardly implemented by requiring \n\ to be large, but the second condition for the curvature 
is more problematic. At this point we have to recall that we are studying isotropic, in par- 
ticular homogeneous, models which are represented by idealized, distributional states in 
the full quantum theory. Such an idealization can lead to problems because one only has 
access to curvature integrated over space rather than local curvature components. In the 
present context, we may have to face infrared problems in the large volume regime because 
the product of curvature and volume of space may be large even if the local curvature is 
small. For instance, if we have a positive cosmological constant A, it will enter the wave 
function in the dimensionless combination Ap which diverges for p — > oo, even though the 
local curvature scale given by A may be small. Similarly, in the positively curved model 
we have the connection coefficient c = \ + jk where r = | comes from the spin connection 
and has the meaning of the integrated intrinsic curvature of space. Therefore, even if the 
extrinsic curvature k is small enough, c may not be so. We will evade those infrared prob- 
lems by assuming c to be small when studying the semiclassical limit. For the flat model 
this can always be achieved by choosing not too large p. 

Furthermore, contrary to a classical symmetry reduction in which a homogeneous ge- 
ometry can be slightly perturbed by adding small non-homogeneous modes, homogeneous 
quantum states are distributional and can only be approximated in the weak topology of 
the kinematical Hilbert space. A consequence is the level splitting in the volume spec- 
trum if we break a symmetry: the simple isotropic volume spectrum becomes increasingly 
complicated as in the homogeneous case or in the full theory. In particular, whereas the 
isotropic volume spectrum has an increasing level distance for large j, the full volume 
spectrum becomes almost continuous (similarly as discussed for the area spectrum in ]16[ ) . 
Such an almost continuous spectrum makes the transition to a classical geometry with 
its continuous volume obvious; but for this also a spectrum with decreasing relative level 
distance is sufficient, as is the case in isotropic models (compare with the equidistant en- 
ergy spectrum of the harmonic oscillator which does not prevent the correspondence to a 
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continuous energy spectrum in the classical regime for large energies). So compared to a 
given volume, the change caused by increasing the time n to n + 1 is always negligible and 
cannot be detected by a classical observer. 

We incorporate this observation in our main condition (pre-classicality f3~ll) for a semi- 



classical regime: the wave function s n ((j)) must not depend strongly on n in the large 
volume regime in order to be regarded as being semiclassical there. More precisely, we 
demand that it is possible to interpolate between the discrete labels n and define a wave 
function ^(a) := s„( a ) with n(a) := 6 sgn(a)a 2 7 _1 Zp 2 with a ranging over a continuous 
range (using \a\ = \/\p\ = y/llp^/\ n \/Q for large \n\ as interpolation points) which varies 
only on scales much larger than the Planck scale. At this point we may have to face the 
above mentioned infrared problems: a cosmological constant leads to a wave function with 
wave length (Aa) _1 which inevitably becomes smaller than the Planck length for large a. 
Generically, there should be a regime in which curvatures are small and the volume is not 
too large in order to allow a continuous time approximation i/j(a). 

Given a wave function -0(a) interpolating the discrete function s n , we can approximate 
the action of the Hamiltonian constraint by derivative operators. The basic operators are 
the difference operator A = 2i sin |c and the mean operator fi = cos |c which have the 
leading order action 

7Z 2 dip 

(As) n ( ) = s n(a )+i - s n(a) _i = 6 P + 0(/ P /a 5 ) (30) 

and 

Mn(„) = l(s n{a )+i + s„( a )-i) = ^(a) + 0(Zp/a 4 ) (31) 

following from a Taylor expansion. The higher order corrections also contain higher deriva- 
tives, but in the semiclassical regime we only need the leading order resulting in the stan- 
dard Wheeler-DeWitt operator 



kH^ ~ -96(*A/2) 2 • fa ~ -6 f-fi 7 Z P JL-^j 



2 

a 



for large a where we inserted A and \i ~ 1 in (2ip, used (|TT|) and expanded the volume 
eigenvalues in j ~ 3a 2 7 _1 /p . This is exactly what one obtains from the classical constraint 

K H^ = -6c 2 sgn(p)v^ 

in standard quantum cosmology by quantizing 3c = — ijlpd/dp. In our framework, how- 
ever, this is only an approximate equation valid for large scale factors where a continuous 
time approximation of the discrete time wave function is possible. Analogously, one can 
show that the term P in the Lorentzian constraint has the correct behavior semiclassically 
which is also true for the spatially positively curved model. For the Euclidean part of the 
constraint, this can most easily be seen in the connection representation where the above 
expansion of the difference operator corresponds to an expansion of sin |c in c. Thus, by 



construction of the Hamiltonian constraint G2J an expansion will result in the standard 
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quantum cosmology expression at leading order. This observation also shows why small c 
are important in the semiclassical regime. 

Arrived at the standard quantum cosmology framework, one can use WKB-techniques 
in order to derive the correct semiclassical behavior. This demonstrates, ignoring possible 
infrared problems, that isotropic models with the evolution equation (|29|) have the correct 



semiclassical behavior for large volume which is achieved in a two-step procedure [36]: First, 



the discrete time behavior has to be approximated by introducing a continuous time and 
interpolating the wave function which results in standard quantum cosmology. In a second 
step, one can perform the classical limit in order to arrive at the classical description. Since 
at fixed k and h the parameter 7 determines the scale of the discreteness, one can also 
describe this by a two- fold limit 7 — > 0, n — > 00 followed by h — > 0. 

In this picture, a universe is fundamentally described quantum geometrically in a dis- 
crete time, and standard quantum cosmology only arises as an approximation which is not 
valid close to the singularity. This explains the large discrepancies of standard quantum 
cosmology and loop quantum cosmology regarding the fate of the singularity. Also the 
issue of choosing boundary conditions, which are usually imposed at a = in standard 
quantum cosmology, appears in a different light. In fact, our discrete time evolution can be 
seen to lead to dynamical initial conditions which are derived from the evolution equation 
and not imposed independently |31[ 



Moreover, using higher order corrections in (|30|) and (|3l|) and in the expansion of the vol- 
ume eigenvalues one can derive perturbative corrections to the standard Wheeler-DeWitt 
operator. Thereby, one obtains an effective Hamiltonian containing higher curvature and 
higher derivative terms. The closer one comes to the classical singularity, the more of 
those perturbative corrections are necessary, until such a description completely breaks 
down at the classical singularity. Because we know the non-perturbative equation which 
is discrete in time, we can see that a perturbative formulation cannot suffice: even if one 
knew all perturbative corrections, it would be very hard to see how they add up to the 
discrete time behavior without knowing the non-perturbative formulation. Furthermore, 
we see that a non-locality in time caused by the discreteness is responsible for higher order 
corrections. This also gives an indication as to why general relativity is perturbatively non- 
renormalizable: adding local counterterms to a local action can never result in a non-local 
behavior like that observed above. 



5 Quantum Flat Space 

As we have seen, the complete Lorentzian constraint is of an awkward form which makes it 
complicated to find explicit solutions. Nevertheless, for the spatially flat model it is quite 
simple at large volume, where the coefficients containing differences of volume eigenvalues 
are nearly identical. The Euclidean part is then of the form of a squared difference opera- 
tor, whereas the Lorentzian constraint is effectively of fourth order. For a complete solution 
we also need to take into account the behavior at small volume, and we will demonstrate 
for the Euclidean part that its values at the classical singularity are crucial for the cor- 
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respondence with the classical situation. Recall that classically we have two solutions in 
the vacuum case (Appendix [A]), a degenerate one given by p = and Euclidean four-space 
characterized by c = (vanishing extrinsic curvature). Standard quantum cosmology leads 
to two independent non-degenerate solutions £i(c) = 5(c) and 62(c) = 5'(c) in the connec- 
tion representation only one of which corresponds to Euclidean space. Moreover, there 
are problems in standard quantum cosmology because it is not possible to invert triad 
operators, even though there is no classical singularity in flat space. 

5.1 Quantum Euclidean Space 

We have already shown that well-defined quantizations of inverse triad operators do exist 
in loop quantum cosmology, so that we do not have to deal with the second problem. The 
first problem of too many solutions will now be investigated for the simplest model, the 
vacuum Euclidean constraint with flat spatial slices. At first sight, it seems to be more 
severe because the discrete Euclidean evolution equation is of order eight, so we have to 
expect eight independent solutions. We already know that one of them is the degenerate 
solution s n = s 5 n fi corresponding to p — 0, and we have to study the remaining seven 
solutions which all lie in the continuous part of the spectrum of the constraint. Since we 
are only interested in solutions with classical regimes we also impose our condition for 
pre-classical behavior, namely that the wave function s n does not vary strongly from n to 
n + 1 for large \n\. This condition, which is the only one besides the evolution equation, 
ensures the possibility of semiclassical behavior for large volume and is independent of the 
explicit form of classical solutions. 
We have the difference equation 

= sgn(n + 4) (Vi| n+4 | - ^| n + 4 |-i) *n+4 - 2 sgn(n) (Vi w - Vi \ n] _^ s n 

+ sgn(n - 4) (Vi| n _ 4| - Vi, n _ 4hl ) s n -4 (32) 

which can immediately be seen to split into independent equations for the four sequences 
Sim, S4m+ij ■54 m+ 2, and S4 m+ 3 with m G Z. The first sequence contains the classical sin- 
gularity at m = and is subject to the consistency condition arising from the vanishing 
highest order coefficient for n = —4. Therefore, it has only one independent solution (the 
would-be other one is the degenerate solution), whereas the three other sequences all have 
two independent solutions because they are subject to a difference equation of order two 
with never vanishing coefficients. By sticking together arbitrary solutions for all four series 
we get all the seven independent non-degenerate solutions of ([32]). At this point, we can 
already impose our selection criterion of pre-classical behavior for large volume: In this 
regime, all four series have solutions S4 m +j ~ aj(4m + i) + 6j, i = ... 3, with an addi- 
tional consistency condition relating ao and 60 which can only be determined in the small 
volume regime. These solutions can be seen by noting that for large \n\ the coefficients of 
s n +4, s n , and s n _ 4 in ( p2[ ) are nearly identical. Our classicality condition then tells us that 
a o — a i — a 2 — &3 =: a because otherwise the complete solution would vary strongly for 
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large \n\ when jumping between the four series: e.g., S4 m+ i — s^ m = 4(ai — ao)m + b\ — b 
becomes arbitrarily large for large \m\ if a\ — 7^ 0. Furthermore, the coefficients 6, cannot 
be very different from each other for the same reason where the difference bi — bj affects 
the amplitude but not the wave length of the variation in the wave function. According 
to our condition that there must not be a variation at the Planck scale, regardless of the 
amplitude, we also have to set bo = bi = b 2 = b 3 =: b (one may permit differences in the bi 
bounded by some small parameter e, but this will only lead to the solution with identical 
b, which is interpreted as the semiclassical part, together with small non-pre-classical con- 
tributions). Therefore, the classicality condition reduces the eight parameters dj, bi to only 
two parameters a and b determining the form s n = an + b for large \n\ without reference 
to a particular classical solution. Now we are in a position similar to that of standard loop 
quantum cosmology: we have two independent solutions only one of which can correspond 
to classical Euclidean four-space. But we still have the consistency condition relating a 
and b which reduces the two independent solutions to only one. To find this unique non- 
degenerate solution with pre-classical behavior we have to use the full equation ([32]) also 
in the small volume regime and in fact right at the classical singularity. This regime is, 
as demonstrated in the preceding section, not accessible to standard quantum cosmology 
which explains the fact that there are too many solutions in this approach. However, it is 
not guaranteed that the unique solution corresponds to the classical solution c = which 
can only be decided when we know its explicit form. Since we need the evolution equation 
at the singularity for a complete solution which also affects the large volume behavior by 
fixing the relation between a and b, we can perform a crucial test of loop quantum cosmol- 
ogy for very strong fields by comparing its unique solution with the classical solution at 
large volume. 

The consistency condition only appears for the sequence s^ m on which we can focus 
from now on; the remaining three sequences are then fixed by the pre-classicality condition. 
We first look at the branch for m > where we choose the only free parameter s^. whereas 
the coefficient so = is fixed (or drops out if non-zero), all other coefficients are determined 
by (^2|) which can be solved for the highest order component 
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For negative n we can do the same, leading to s_4, s_8, ... in terms of S4. The result 
is (p3|), now with arbitrary m G Z, which is the exact solution for the sequence Si m and 
automatically fulfills the consistency condition since we started from so = 0. Implicitly, 
(j33]) determines the relation between a and b such that we now also have the unique solution 
with pre-classical behavior given by 



V 2 -V x 



l M y " «*■ (34) 



1 



Dropping constant factors (or choosing S4 := 2(X^ — Vi) 1 ), this yields in the connection 
representation 

00 9-1-1 

lf){c) = =X) y.;_y. 0(c) (35) 

where the expression 2j + l = \n\ (rather than another linear function in j) appears because 
the consistency condition fixed the relation between a and b. This allows us to check 
the compatibility of this condition, which arose because of the structure at the classical 
singularity, with the correct classical behavior: we expect a solution which is related to the 
5-function in c incorporating the classical solution c = (see also App. |A]). 

More precisely, since we chose an ordering with the triad components to the right in 
the connection representation, which has been seen to be necessary in order to remove 
the singularity, we should expect a quantization of a = sgn(p)y]p[ to map if) to the 
delta function (the classical constraint is — 6c 2 sgn(p)y / [p| in the ordering chosen for the 
quantization). A possible quantization of a maps 

a( j = 2i( 1 l 2 P )- 1 (V j+l2 -V j _i) Xj 

with an analogous action on \ji which has the correct large-j behavior for \/\p\ and which 
maps Q to ixj (and to zero for j = — |). The last property is necessary due to the sgn(p) 



which can be seen from the action on |n). One can also use the techniques of |25[ in order 
to derive this quantization: writing 

sgn(p)/H = a = |Af4 = -|tr(A 7 a}r,) 
and using Thiemann's quantization of the co-triad components we have 

a = -|i( 7 / 2 )- 1 ^tr(A / /i / [/i7 1 , V\) = -4i( 7 / 2 ,)- 1 (singe) V cos(i C ) - cos(fc) V singe)) 
1 

acting as 

a Xj = -2<( 7 ^)- 1 (V i+ i - ^._i)G 

aQ = 2i( 1 l 2 p y 1 (V jH -V j i) Xj . (36) 



The eigenvalues of a are 



a 



, - ±2{ 1 ll)-\V 3+l -V 1 _ h 
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which for large j is 

Applying a of fl36|) to our solution (|35|) yields 

a^(c) = 2z( 7 ^)- 1 £(2j + l) Xj {c) = 2i( 1 ll)- l 5(c) 

3 

which in fact is proportional to the 5-function on the configuration space 577(2). For 
this it is crucial that we have the factor 2j + 1 which appears uniquely only if one uses 
the consistency condition arising from the behavior at the classical singularity. Thus, 
the unique solution to the Euclidean constraint with semiclassical behavior in spatially 
flat isotropic loop quantum cosmology correctly corresponds to Euclidean four-space with 
vanishing extrinsic curvature. Similarly in other models, the consistency condition together 
with the pre-classicality condition on the variation of the wave function for large a always 
selects a unique solution. In this way dynamical initial conditions are derived from the 
evolution equation |y] and not imposed ad hoc as usually done in standard quantum 
cosmology. 

5.2 Lorentzian Constraint at Large Volume 



For large \n\ the Lorentzian constraint equation (|29f) simplifies since V\ n+ k\/2 — V(| n +fc|_i)/2 
is approximately independent of k for n 3> k and the extrinsic curvature coefficients k^ 
are nearly one. Introducing t m := S2 m , the vacuum constraint equation takes the form 

= |(1 + 7" 2 )W4- W2 + |(3-7 _2 )im-t m -2 + |(l + 7" 2 )i m -4 

= i[(A 4 + 4 7 - 2 AV)i 

= i[A 2 (A + 2z 7 -V)(A-2* 7 -V)t] m 

using again the central difference and mean operators A and \i. Because they commute we 
can split this equation into 

A 2 t = 
or (A + 2z7~V)* = 
or (A-2i 7 ~V) i = 

with independent solutions 



= ci t {2) 

t (3)_„ f l-ij-' Y ,(4) _ fl~^- 1 



1 + tj 1 ) \1H-Z7 1 
For 7 = 1 we can also use the form 

= c 3 Re(i m ) = cos(mf ) , t$ = c 4 lm(i m ) = sin(mf ) 
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As with the Euclidean constraint we have still more solutions since the original equation 
( p9|) is of order 16. But again it splits into four sequences with solutions as above which 
can be put together to form all 16 independent solutions. Thanks to the pre-classicality 
condition of mild variation only four of them are relevant. However, only two, t^ 1 ' and 
can correspond to the two solutions of the second order standard Wheeler-DeWitt 
equation; the rest has to be excluded on general grounds which again are provided by 
pre-classicality. This immediately excludes and if 7 = 1 which jump between the 
values and ±1, and also for 7 7^ 1 if 7 is of order one or less (which is true for the physical 
value |37| |38|]): tm and tm are of the form exp(±im9) with 9 = arccos[(l — 7 -2 )/ (1 +7~ 2 )] 
which is 7r > 9 > I for < 7 < 1. 

We now arrived at two independent solutions allowed for a semiclassical analysis from 
which a particular combination is selected by the consistency condition at the classical 
singularity. To evaluate this we would need to take into account the exact equation (f29|) also 
for small volume, from which we refrain here. This appears not to be possible analytically 
in closed form, but can easily be done in a numerical study. Nevertheless, it is easy to see 
that the Hamiltonian constraint equation has solutions of the same semiclassical behavior 
for all values of 7, which contradicts the hope that the Hamiltonian constraint equation 
selects a value for 7. We can only conclude that a large parameter 7 ^> 1 would lead to 
additional pre-classical solutions lacking any correspondence to a classical solution (e.g., 
oscillating solutions for flat space). This may be taken as an argument that 7 cannot be 



much larger than one, in coincidence with [37, 3S 



6 Comparison with Other Approaches 

A resolution or avoidance of the classical singularity has been claimed before in a variety of 
approaches. After the singularity theorems of general relativity had been established, 
it became clear that one has to couple matter which violates energy conditions in order 
to evade them (it is sufficient to violate only the strong energy condition ||). This can 
be achieved with either classical @ or quantum matter || |J] leading to a bounce in the 
evolution of a universe at positive radius. However, this conclusion is model and parameter 
dependent and, therefore, not a generic behavior. Furthermore, the threat of a singularity is 
still present in gravity, but only avoided by particular types of matter (and, since quantum 
matter field theories have their own divergences, it may be dangerous to call upon quantum 
matter for a rescue from the singularity). 

Another idea to evade the singularity theorems consists in changing general relativity. 
Since its action is deemed to be only an effective action of something more fundamental, 
there can be correction terms being non-linear in the Ricci scalar R. Inclusion of the lowest 
order correction quadratic in R has been shown to yield solutions which do not encounter a 
singularity || [7|, ||. Again, the conclusion is parameter dependent and, in fact, inconsistent 
because the non-singular solutions emerge only as artifacts of the truncation of the higher 
order corrections M. This situation suggests that a complete non-perturbative formulation 
is necessary for an investigation of the fate of the classical singularity. 
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The original approach to this problem in canonical quantum gravity, which is non- 



perturbative, was started in [|T0[ . Here, it was proposed that one should use the boundary 
condition t/*(0) = of a vanishing wave function right at the classical singularity It has 
also been speculated that this condition is enforced by an ad hoc Planck potential which is 
relevant only for small scale factors ||40|| . However, as argued in ||12|| , this requirement by it- 
self cannot be regarded as a sufficient condition for the absence of a singularity because the 
scale factor has a continuous spectrum in standard quantum cosmology (note that in loop 
quantum cosmology the spectrum of the scale factor is discrete, but the mechanism which 
removes the classical singularity is very different from De Witt's proposal); in addition an 
appropriate fall-off behavior of i[) close to a = would be necessary. In this context, it has 
also been suggested, sometimes for purely mathematical reasons in order to obtain a self- 
adjoint id/ dp, to extend minisuperspace to include negative volumes, which would remove 
the boundary at a = and allow wave functions to extend into this regime [34]] (although 
this may seem similar to the evolution through the singularity derived in the present pa- 
per, it is not to be confused with our negative p which still leads to positive volume: in 
contrast to negative definite metrics, negative triads are allowed classically even though 
disconnected from the p > sector if one requires non-degeneracy). The negative metric 
branch lacks a classical interpretation, and so the wave function is completely quantum 
without semiclassical interpretation in this large region of the configuration space. It has 
been suggested that the transition to negative volume should be interpreted as a signature 
change to Euclidean space-time or a "tunneling from nothing" |f42| . 

As derived in this paper, loop quantum cosmology is able to describe the behavior of 
a universe close to the classical singularity. It always leads to a decoupling of Sq from 
solutions of cosmological behavior, reminiscent of DeWitt's i/j(0) = 0. However, since 
from the present perspective standard quantum cosmology completely breaks down in 
this regime, a condition for the standard wave function ip{a) is no longer meaningful at 
a = 0. This condition for sq is in fact important for the absence of a singularity, and it 
serves to select a unique superposition of the two WKB components when evolved into the 



semiclassical regime pi] . At a = we have a transition to another branch which opens 
to large positive volume at negative time (rather than a "tunneling from nothing", our 
picture of an evolving universe could be described as "tunneling through nothing" if one 
wants to identify "nothing" with the degenerate state |0)). 

This picture of a branch preceding the classical singularity is reminiscent of the "pre- 
big-bang" scenario (and other, more recent constructions) of string cosmology [44 



where a contracting universe preceding the singularity has been claimed which should be 
connected to the present expanding branch through a high curvature regime. However, 
lacking a non-perturbative framework, this claim cannot be substantiated at present. 

The possibility of an evolution through a degenerate state also reminds of results which 
have been obtained in the context of mirror symmetry: by mapping the degenerate state to 
an equivalent non-degenerate one it is possible to extend the evolution through a singular 
geometry E5|, However, this has been demonstrated only for very special spaces 

where one generally focuses on the geometry of compact directions; these spaces are not 
sufficiently realistic to be of direct physical interest, let alone cosmological models. 
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7 Conclusions 



A reduction of quantum geometry to isotropic geometries leads to models in which explicit 
computations are possible, thanks primarily to the completely known volume spectrum. 
Therefore, they provide an ideal test arena for the techniques of quantum geometry, which 
turn out to work without any problems. Moreover, they are interesting in their own right 
as cosmological models where they shed light on certain aspects of the classical singularity 
which were not illuminated in any other approach to (quantum) cosmology. A kinematical 
indication for a better behavior of loop quantum cosmology close to the classical singu- 
larity comes from a quantization of the inverse scale factor p5[ using techniques for the 
quantization of matter Hamiltonians to densely defined operators in quantum geometry 
|2~7fl . The result is a bounded inverse scale factor which does not diverge even if the volume 



is zero, a result which is possible only in a quantum theory of geometry. In the present 
paper we established the absence of a singularity by studying evolution equations at the 
dynamical level: whereas classically the evolution breaks down at the singularity, in loop 
quantum cosmology we can evolve through it. Since evolution through a degenerate state 
is possible, one could also obtain topology change in quantum geometry. 

At this point we explain in more detail in which sense the singularity is absent in 
loop quantum cosmology. One might think that there is still a singular space geometry 
of vanishing volume. But this is not as problematic as in the classical theory since, e.g., 
the inverse scale factor does not diverge. One should also keep in mind that vanishing of 
volume is possible even classically without a singularity: it may just signal the presence of 
a horizon as is the case in (non-singular) de Sitter space-time when sliced by flat spaces. 
In this case one can, of course, evolve through the horizon by choosing an appropriate 
time coordinate. In contrast to a singular space-time, such a manifold is not geodesically 
incomplete. On the other hand, the existence of an incomplete curve together with energy 
conditions for the matter inevitably leads to a curvature singularity At such a point, the 
curvature tensor cannot even be interpreted in a distributional sense and so Einstein's field 
equations break down. There is then no means to extend the singular space-time beyond 
the singularity in a unique manner (it may be possible to extend a space-time continuously, 
but never uniquely). The evolution equation (|2"5| ) of isotropic loop quantum cosmology, on 
the contrary, never breaks down and so always gives rise to a unique extension through the 
quantum regime containing the classical singularity. 

In deriving this behavior it was our strategy always to be as close to the full theory 
as possible, e.g. when quantizing the Hamiltonian constraint. Although one might have 
simplified some expressions in a model-dependent way, this would have lead to deviations 
from the methods of the full theory. For instance in the case of flat spatial slices the 
classical Lorentzian constraint and its Euclidean part only differ by a factor 7 -2 so that 
one might be tempted to use this relation also in a quantization which would strongly 
simplify the quantum constraint. However, this fact crucially depends on i) K l a oc A\ 
(due to = 0) and ii) 6ij k A % a A{ oc F^ b (due to dA = 0) both of which fail in other 
homogeneous models, let alone in the full theory. Making use of these relations would 
simplify the computation, but the results were not trustworthy since the contact to the 
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full theory would have been lost. Note that in standard quantum cosmology one does not 
have a corresponding quantization of the full theory as guidance, and thus lacks means to 
evaluate manipulations. In fact, loop and standard quantum cosmology differ from each 
other right in the regime where quantum gravitational effects are important. Standard 
quantum cosmology uses quantum mechanical methods and phenomena (like the tunneling 
effect) but not full quantum gravity, whereas loop quantum gravity is close to the full theory 
of quantum geometry. Some of its quantization techniques may seem to be unfamiliar from 
the viewpoint of quantum mechanics, but they are necessary since analogous techniques are 
required for a consistent quantization of gravity. Moreover, quantization ambiguities (like 
the one described at the beginning of this paragraph) are severely restricted by requiring 
that analogous quantizations must be possible in the full theory. 

For the absence of a singularity the form of the matter coupled to gravity is irrelevant 
because the removal of the singularity is completely due to quantum geometry. On the 
other hand, the (non-symmetric) factor ordering of the Hamiltonian constraint is crucial 
for this result, so that demanding a non-singular evolution in quantum cosmology fixes the 
factor ordering ambiguity of the constraint: the scenario derived in this paper is possible 
with only one ordering, which belongs to the three standard choices (one could still choose 
different orderings of the extrinsic curvature operator entering the Lorentzian constraint, 
but since it is an observable in the kinematical sector it should be ordered symmetrically 
as done here). 

Close to the classical singularity the discrete structure of space and time in quantum 
geometry is important which leads to large deviations from standard quantum cosmology 
[TIJ. This framework arises here as an approximation which is good only at large volume 



where the discrete volume spectrum is washed out to a continuous spectrum by inaccura- 
cies. However, the exact description of loop quantum cosmology is also necessary to fix 
a unique solution which can be seen in the explicit solution corresponding to Euclidean 
four-space. For cosmological models with matter there is still a unique solution with appro- 
priate semiclassical behavior: initial conditions are not imposed ad hoc but instead derived 



from the dynamical laws [[31]. Taking into account the discreteness of the spectrum one 
can derive perturbative corrections for an effective Hamiltonian of standard quantum cos- 
mology, but the completely non-perturbative description with discrete time is needed in 
order to study the fate of the classical singularity. 

As an intuitive picture of the evolution of a universe in loop quantum cosmology we 
obtain the following one: starting in a semiclassical contracting state, it reaches a degen- 
erate stage seen as the singularity in classical cosmology, in which the universe bounces 
off in order to enter an expanding branch. The further fate, whether it expands forever or 
recollapses in order to start a new such process, depends on the matter content and the 
value of the cosmological constant. 
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Appendices 



A Euclidean Space in Standard Quantum Cosmology 

Here we present classical and standard quantum cosmology of flat Euclidean space. For 
T = the constraint equation (|i~3| ) becomes = — Qk^c 2 sgn(p)^^/\p\ = having two 



solutions p = or c = 0. The first one only appears in a triad formulation and is of 
no physical interest because space is completely degenerate. Using fll5]) one can see that 
the second solution requires the extrinsic curvature to vanish, whereas it leaves the metric 
arbitrary, and thus gives flat Euclidean space (the metric does not change in coordinate 
time because of {p, H^} « on the constraint surface). 

Standard quantum cosmology proceeds as follows: Using a factor ordering in which 
the metric variables appear on the right (this is the ordering resulting from loop quantum 
cosmology) solutions in the c-representation are such that a quantization of sgn(p)y|p| 
acting on the wave function yields either 5(c) or S'(c) due to the factor c 2 . There are some 
problems already in this simple model: First, sgn(p)^^/\p\ cannot be quantized to an invert- 
ible operator when the range of p contains the value zero, and a procedure like that leading 
to a bounded inverse scale factor in loop quantum cosmology is not available in standard 
quantum cosmology. Although Euclidean space does not contain a singularity, the point 
p = leads to problems in the quantization. Second, there are two independent solutions 
only one of which, 6(c), corresponds to the classical solution; and there is no independent 
argument to exclude the other solution without referring explicitly to the classical situa- 
tion. Both problems are solved in loop quantum cosmology, which is intimately related to 
the fate of the classical singularity. 



B Hamiltonian Constraint for Models with Positive 
Spatial Curvature 

For the isotropic, spatially positively curved model we need to take into account the holon- 
omy hu^j\ in flUf ) which is Kk if &ijk is positive and if ejjk is negative. With this we 
have for isotropic holonomies 

Y,^ JK ^{h I hjh^ 1 hfh^ J] h K [h K \V]) = 3tT^h 1 h 2 h^ 1 h^ 1 [h^ 1 ,V}^ 

UK 

-3tr ^hxh^h^hl^ 1 ,V]\ 

which yields 

Hf> = -48i( 7 ^ 2 )- 1 (sin(§c) - 2 sin 5 (i C ) - 2 sin 2 (i C ) cos 3 (i C )) 
x ( sin(|c)Ucos(|c) — cos(|c)Usin(|c) j 
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with action 

Jff>|n) = 3(7^)- 1 sgn(n)(yi w -yi w _ 1 )(i(l + i)|n + 5) + |(l-5i)|ri + 3) 
- (1 - i)|ra + 1) - (1 + i)|n - 1) + |(1 + 5i)\n - 3) + ±(1 - i)|ra - 5)) 



From this we obtain the extrinsic curvature operator 

5 



K + \n) = \ll £ ^f„ 9) |n + 9> 

q=-5-q odd 
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= ±6(5 T 0( 7 /|)- 3 (V |W 


- Vi H 


-l) (^(1^31- 


-1) 


- v m 


-1) 


^(±5) 
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This leads to 

5 

^sin(|c)i^ + cos(|c) — cos(|c)i^ + sin(|c) j |n) = |i/p £ k^^\n + q) 

q=-5;q odd 

with 

1.(9) _ l(lf(<l) _ Jf(<l) \ 
h +,n ~ 2\ +,n+l JX +,n-l) 

which is non-zero for all n and q (K^ n is zero if and only if n — 0). 
Taken together, this yields 

P + \n) = -Ul + l-^llr^gnin) (v^-V^) £ A[;^\n + 2k) 



k=-5 

with 

4(0 ■= V fc(< , 



<j+r=Z ; -5<5,r<5 ; <j,r odd 
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